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Abstract 

In this paper we determine a formula for calculating the refractive index n for the acoustic 
equation from the partial Dirichlet to Neumann map(DN) associated to n. We apply these 
results to identify locations and values of small volume perturbations of this refractive index at 
fixed frequency u. 

Identification de I'indice de refraction de I'equation acoustique a frequence fixe 

Resume 

Dans cette Note nous determinons une formule pour calculer I'indice de refraction n pour 
I'equation acoustique a partir de I'application Dirichlct-to-Neumann partielle (DN) associee 
a n. Nous nous appliquerons ces resultats pour identifier les locations et les valeurs des petites 
perturbations volumiques associees a cet indice de refraction pour une frequence to fixe. 

Version frangaise abregee 

Cette Note traite un probleme inverse pour I'equation acoustique avec frequence fixe. Le but est 
d'identificr I'indice de refraction associe apres avoir deriver une formule appropriee a I'aide de 
I'application Dirichlet-to-Neumann partielle (DN) associee a n et de donner plus explicitement 
cette identification. Concernant I'utilisation de I'application Dirichlet-to-Neumann partielle en 
probleme d'identification, nous pouvons trouver, par exemple, le travail De Kohn et Vogelius [5]. 

Soit fl C R"^ un domaine borne avec un bord de classe C^. Nous noterons v la normale 
unitaire sortante du bord dfl. Soit T une partie ouverte reguliere de la frontiere dfl. Supposons 
que fl contient un nombre fini d'inhomogenetes, chacune de la forme zj + aBj, ou Bj C R^ 
bornee et contenant I'origine. 

Soit n(a;) € C°(ri) I'indice de refraction non perturbc. Nous supposons que n(x) est connue 
sur un voisinage de dft. Notons par nj{x) G C^{zj + aBj) I'indice de refraction de la j-eme 
inhomoe'ete, Zj + aBj. 

Considerons I'equation acoustique, a frequence fixe, en presence d'inhomogeneetes 

(A -I- uj^na)ua = et n 

et definissons I'application Dirichlet-to-Neumann partielle associee a Ua par: An^ (/) — 

du 

-^-^|r pour tous / G H^iY), Ici _ff 2 (r) designe I'espace trace. 

Bukhgeim et Uhlmann 4^ ont utilise les solutions complexes de I'optique geometrique pour 
prouver que la connaissance des donnees partielles de Cauchy pour I'equation de Schrodinger 
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determine le potentiel de maniere unique sachant que le potentiel q e L°°(S7) est connue dans 
un voisinage du bord. Nous pouvons utiliser leurs travaux pour montrer notre version physique 
et nous utihsons ainsi ces solutions complexes de I'optique geometrique construites dans leur 
travaux pour prouver notre niethode de reconstruction. 

Ensuite, nous derivons une formule de calcul de I'indice de refraction n a partir de I'application 
Dirichlet-to-Neumann partielle associe a n et F en modifiant la procedure de reconstruction de 
Nachman Enfin, nous profitons des proprietes des noyaux des operateurs ainsi introduites 
pour presenter plus explicitement cette formule de reconstruction. Nous considcrons le probleme 
inverse d'identification des lieux de petites perturbations volumique de I'indice de refraction et 
nous appliquons nos resultats pour reduire ce probleme inverse au calcul de transformation de 
Fourier inverse. 



1 Problem formulation 

The paper deals with an inverse problem for the acoustic equation with fixed frequency. The 
object is to identify the associated refractive index after deriving a suitable formula by using 
the partial Dirichlet to Neumann map(DN) associated to n and to give more explicitly this 
identification. Concerning the use of the partial Dirichlet to Neumann map in recovering problem 
we can find, for example, the work of Kohn and Vogelius [3- 

Let ri C R'^ be a bounded domain with boundary. We denote by v the unit-outer normal 
to dil. Let F be a smooth open subset of the boundary dfl and F^, denotes dil \ F. We suppose 
throughout that R'^yfi. Introduce the trace space 



H^T) ={/ e H^dn)J = on Tc]. 



Here and in the sequel we identify / defined only on F with its extension by to all dfl. It is 
known that the dual of i?^(F) is ^^-^(F). 

Assume that n contains a finite number of inhomogeneities, each of the form zj + aBj, where 
Bj C R^ is a bounded, smooth domain containing the origin. The total collection of inhomo- 
geneities is Ba = ^^^li^j + ctBj). The points zj S fi, j = 1, . . . , m, which determine the location 
of the inhomogeneities, are assumed to satisfy the following inequalities: 

(1) \zj - zi\> Co > 0,y j ^ I and dist(zj,afi) > Co > 0,V j, 

where Cq is a positive constant. Assume that a > 0, the common order of magnitude of the 
diameters of the inhomogeneities, is sufficiently small, that these inhomogeneities are disjoint, 
and that their distance to R'^ \ is larger than cq. 

Let n(a::) g C°(il) denote the unperturbed refractive index. We assume that n{x) is known on 
a neighborhood of the boundary di^. Let rij{x) € C^{zj + aBj) denote the refractive index of 
the j-th inhomogeneity, Zj aBj . Introduce the perturbed potential 



(2) na(a;) 



n(a;), a; e f2 \ Sq, 

x\.j{x), X (z Zj + aBj, i = 1 



Consider the acoustic equation, at fixed frequency, in the presence of the inhomogeneities 

(A + uj'^na)ua = in 

uc.\an = f eSHr), 

and define the local Dirichlct to Neumann map associated to Hq, by :An„(/) = ^^|r for all 
/ e H2{T). Where 775 (r) is the trace space 



Bukhgeim and Uhlmann [4J used the complex geometrical optics solutions vanishing on the 
complementary of F to prove that the knowledge of the partial Cauchy data for the Schrodinger 
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equation on any open subset T of the boundary determines uniquely the potential q provided 
that q S L°^ (ri) is known in a neighborhood of the boundary. We may refer to their work to 
show our physical version and we use these complex geometrical optics solutions, constructed 
in their work, to prove our reconstruction procedure. 

Next, we derive a formula for calculating the refractive index n from the partial Dirichlet- 
to-Neumann map (DN) associated to n and F by modifying the Nachman's reconstruction 
procedure [5] . It turns out that if the refractive index is a-priori known in a neighborhood of the 
boundary the derivation of a reconstruction formula are surprisingly simple. Finally, we profit to 
properties of the kernel of the introduced operator to explicit more this reconstruction formula. 
We consider the inverse problem of identifying locations of small volume fraction perturbations 
of the refractive index and we apply our results to reduce this inverse problem to calculations 
of inverse Fourier transforms. 



2 Identification procedure 

Let no e L°°{il) be a known function. Assume that n = no almost everywhere in a neighborhood 
of dfl. We extend n and no by in R'^. Let p e R'^ \ {0} and define Up to be the solution of 

Au„ = in R^ \ n, 

(3) 

(^A + n)up = in n, 
we find that Up|r solves the (hyper-singular) integral equation on the open surface F: 

(4) AnK|r)+/ ^ (x,y)uMy)ds{y)^p-^{x)e^-'', V x e F, 

/r dv{x)dv{y) 

where is a well defined exterior Dirichlet Green's function for A. 

Now, according to [1, and 3J we can prove the following global uniqueness result associated 
to refractive index in absence of any inhomogeneities (in presence of the inhomogeneities we 
may obtain a similar result), as follows. 

Proposition 1 Suppose that a = 0. Let n^ real-valued in L°°{Q),i = 1,2. Assume ni = n2 
almost everywhere in a neighborhood of the boundary dfl and Am = An^ ■ Then ni = n2 almost 
everywhere in O. 

Proof. Let —1<S<0 and introduce the weighted L^-space 

LliR') = {/ £ LUr')\ l^^il + \x\y\f{x)f dx < +^}. 
From [7] , we know that if n^ is an extant function of n^ defined by: 



ni{x) 



ni{x), X efl, _ 

0, xe R^\n, 



then the solutions of (^A + hi)vi = on R'^ can be given by: 

(5) v^ = e''■P^il + ^n,ix,p^)),^ = l,2 

for \pi\ sufficiently large with ipn-{-,pi) € L^(R^). Moreover, we have 

(6) ll^n,(-,P.)IL.(R3,<^. 

Now, we set 



f] .,k + 1. , 77 .,k ~L 

Pi^2+^{^^) and p2 = +H^— ) 
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where r],k,l e such that rj-k^rj-l^k-l^O, and lyyp = |fcp + l^p. 
Let fi' CC fl, fl' open set with boundary and n\n' is connected. Define 

N{n) e I (^A + n)i; = in r2, w = on r^j 

and 

Af(f7) =|w e H^{n) I (-^A + n)u = in f^j. 

Then, according to [3] the set A^(il) is dense, in the L^{il') norm, in N{il). On the other hand, 
by Green's theorem we have 

(7) / (ni - n2)uiM2 = /(^^W2 - ui^^) lis(x), 

where ds{x) denotes surface measure and ui,U2 G NiVl). Let zi G H^iVi) n N{n) such that 
U2\t = zi\r- The hypothesis Am — gives that 

(8) 2i|r, = 0, zi|r = W2|r ^ ^|r = ^|r- 
Combining relation ([5]) with ([7]), we deduce: 

/ {ni - ■n2)uiU2 dx = [ {^^U2 ~ ui^^) ds 
Jn' Jr di' dv 



dui dzi 
(— — zi - Mi^j ds. 



By Green's theorem, we have: 



[ i^^r-zi - ui^^) ds ^ [ (ni - ni)wizi (ix = 0, 



which imphes 

(ill — n2)uiU2 dx = 0. 



/ 



Now, using density's argument we can approximate any Zi G N{^) by elements of N{^). There- 
fore the functions vi defined in ^ and belong to N{^) satisfy 

/ (ni — n2)ziZ2 dx — 0. 

Letting |Z| — > +cx), and using the estimation ([6]), we deduce: 

(ni'^n2)(/s) = VfcGR'^ 

which achieves the proof. □ 
Next, define the double layer potential 

^pif^^f n f?/. A x,y)f\r{y)dsiy), v/Gi/^(r), 

Jr dv{x)dv{y) 

and set: 

■q ..k + 1. 

where r],k,l £ such that rj ■ k — rj ■ I = k ■ I = 0, and |?7p = |fcp + The following holds. 
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Lemma 1 Assume that is not a Dirichlet eigenvalue of (^A + n) in ft and pi be given as 
in 0). Then, there is a unique solution Wpjr € Hi{T) of ^ such that 



UpAr = (A„ + iVpJ-i(pi • iyix)e^-PAr). 



Next, we can prove the following reconstruction formula by using Proposition [TJ 

Proposition 2 Let no € L°°{fl) be a given function. Assume that is not a Dirichlet eigen- 
value o/(^A + n) in and n = no almost everywhere in a neighborhood of d^. Then 

{n^n^){-k) = \ lim /(A„ + 7VpJ-i(pi • Kx)e^-''Hr)(A„- 



|;|^+oo ' 

A„„)(A„„ +7VpJ'i(p2 •Ka^)e^-''^|r) ds{x), 
where k, pi and p2 are given as in 

Proof Let pi e \ {0} [i = 1, 2). Let Upjr e ^^(r) be the solution of (g]) and define 

On{x,pi) = e'^'''P'Up^{x) - 1, 
then we have Up^ = e~^'''^{l + On{x,pi)). 

We decompose ps = + irj) {i^ = — f ) with ^, 77 S R"^, |^| = |?7| = 1. Then we can write: 

Ap0„ = nxme^ in R3\r„ 

where x(^) is the characteristic function of Q.. In a small neighborhood V{d^) of the boundary 
9f2 let denote the normal coordinate and (/^(j^) a smooth cut-off function which vanishes on 
V{dVt). The fact that fi' is compactly supported in 17 as mentioned in the proof of Proposition 
[11 we have iy9(j^)6'n = 0^ on il' for \p\ sufficiently large. In order to get the aid of Proposition 
[T]we, firstly, demonstrate the following estimate 

Il^n(-,P)||L^(0') < 

To do this, we may set 9^ = <(5(j^)^n in R^ then 

Ap^n = nx{^)h + • • + eJj^ + A(^ + 2-^ • V^) in R^ 

\p\ \pr \p\ 

Obviously that 9n E Lg, therefore by classical results [TJ the following estimate holds 

(10) ||^~n||^.(R3)<^, 

for some positive constant C independent of t. Next, as done in ([9]) we set 

7] .,k-l 

where ri,k,l e R^ such that r]-k~rj-l = k- l = 0, and = |fcp + \l\'^. 
By applying last results, we construct 

Up^ e''-P'{l + enix,pi)) inn' 

and 



v = Vp = e--PHl + en,ix,p2)) 
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to obtain that 

(11) (n-no)(-A;)= lim u p{An ~ Ano)v p ds. 

Finally, by Lemma[T]one can see that the boundary values of the solutions Up|r can be recovered 
from An-The result is then proven. □ 
The following lemma is useful to give more explanation to our reconstruction method. 

Lemma 2 Let pi and p2 be given as in the following asymptotic behavior holds : 

1 



Np^ = \l\L, + 0{-^) as |;| ^+oo, 



where, for i G {1, 2}, is a well defined integral operator on H ^ (F). 

Proof. We give a brief proof for i = I, and the case i — 2 can be given similarly. Recall that the 

kernel of the operator iVp, is J'\ {x, y) where g^ (x, y) = e'^'f^G^ {x, y) and {x, y) is 

av^xjaviyy) ^ ' ' 

a solution of the following integral equation of the first kind : 

f dG^ 

(12) -Gp,{x)^ Gp,{x,y)^^Ax,y)ds{y). 



on 



dv{y) 



Here, the function Gp^ is given by 

On the other hand, one can use dH) to find that 



Pi-C = lC|cos(pi,0 
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1^1 



Therefore, to continue with the proof we may insert the last formula into equation (jl2p and we 
follow the convenable procedure. □ 
Applying Lemma [2] to the results found in Proposition [2l we can prove the following recon- 
struction formula. 

Theorem 1 Let Uq G L°°(r2) be a given function. Assume that is not a Dirichlet eigenvalue 
of (;^A + n) in Q and n = no almost everywhere in a neighborhood of dVL. Then 

(n^no)(-fc) = -4 / ^(^•K^)|r)'e*^ ds(x), 

where r],k £ R'^ such that rj ■ k = 0. 

Now, we apply Proposition [2] and the asymptotic result given in [2] for identifying efficiently 
the locations {zj}JLi of the small inhomogeneities Ba from the knowledge of the difference 
between the local DN maps A„^ — A„ on F. The following result holds. 

Theorem 2 Suppose that we have (QJ), let n G C'^{il), be given by (0). Assume that is not 
a Dirichlet eigenvalue o/ A + n^) in f2 and Ua = n almost everywhere in a neighborhood of 
dfl (for a sufficiently small). Let k, I, pi and p2 be as in Then, the following identification 
holds: 

{n:^n){-k)^\ lim ( {K^ + Np,)-\p^ ■ v{x)e''-P'\r){K^- 

r> 771 



An)(A„ + Np,Y\p2 ■ v{x)e--P-^v) ds{x) = — ^(n(z,) - n,(z,))|S,|e*'-^^- + o{a' 
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As done in Theorem [l] we can deduce the foUowing more explicit result. 

Corollary 1 Suppose that we have the hypothesis of Theorem [H Then, the following identifi- 
cation holds: 

m 

Jt \v\ ^.^1 

where r],k € R'^ such that rj ■ k — 0. 

By neglecting the remainders o(a'^) in Corollary [1] the locations {zjjJLi are obtained as 
supports of the inverse Fourier transform of 

-V2 / ^(r,.Kx)|r)V'=- ds{x). 
Jr \v\ 

If, we get the points {zjYfLi, the values {nj{zj)}jLi could be obtained by solving a linear system 
arising from Corollary [TJ 
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